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Abstract Utilizing a BzK-selected technique, we obtain 14550 star-forming galaxies 
(sBzKs) and 1763 passive galaxies (pBzKs) at 2 ; ^ 2 from the /C-selected (Xab < 22.5) 
catalog in the COSMOS/UltraVISTA field. The differential number counts of sBzKs and 
pBzKs are consistent with the results from the literature. Compared to the observed results, 
semi-analytic models of galaxy formation and evolution provide too few (many) galaxies at 
high (low)-mass end. Moreover, we find that the star formation rate (SFR) and stellar mass 
of sBzKs follow the relation of main sequence. Based on the HSTAVide Field Camera 3 
(WFC3) F160W imaging, we find a wide range of morphological diversities for sBzKs, from 
diffuse to early-type spiral structures, with relatively high M 20 , large size and low G, while 
pBzKs are elliptical-like compact morphologies with lower M 20 , smaller size and higher G, 
indicating the more concentrated and symmetric spatial extent of stellar population distribu¬ 
tion in pBzKs than sBzKs. Furthermore, the sizes of pBzKs (sBzKs) at z ~ 2 are on average 
two to three (one to two) times smaller than those of local early-type (late-type) galaxies with 
similar stellar mass. Our findings imply that the two classes have different evolution modes 
and mass assembly histories. 
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1 INTRODUCTION 

The formation and evolution of massive galaxies (M* > 10^° Mq) at z 2 is a hot issue of observational 
astronomy. There are many reasons, for instance the population of Hubble sequence galaxies is already in 
place at z 1.5-2 (Fang et al. 2012), the universe star formation rate density (SFRD) peaks at z ~ 2 
(Oesch et al. 2012), the specific SFR (sSFR) evolves weakly at z > 2 (Gonzalez et al. 2014), the galaxy’s 
mass grows quickly at 1 < z < 3 (Ilbert et al. 2013), luminous infrared galaxies (Lg-iooo Atm > 10^^ Lq) 
are more common at redshift z ^ 1-3 (Murphy et al. 2013), and the number density of quasi-stellar objects 
(QSOs) has a peak at z 2 (Richards et al. 2006). 

Within the past decade, many novel techniques have been applied to select a sample of massive galaxies 
at the epoch of z 2 and important investigation has been made in our understanding of high-redshift 
galaxies (Chapman et al. 2003; Franx et al. 2003; Daddi et al. 2004; Kong et al. 2006; Dey et al. 2008; Huang 
et al. 2009; Fang et al. 2012; Wang et al. 2012; Fang et al. 2014). Such as submillimeter galaxies (SMGs 
with F’(850 p,m) > 0.5 mJy), distant red galaxies (DRGs with (J — ifjvega > 2.3), ultraluminous infrared 
galaxies (ULIRGs; Ls-iooo Atm > 10^^ Lq), dusty-obscured galaxies (DOGs with (i? — [24])vega > 24), 
etc. Based on a simple two-color {B — z and z — K) approach, Daddi et al. (2004) introduced the criteria 
of (z — K)ab = 2.5 and BzK = (z — K)ab — (B — z)ab = —0.2 to select the sample of z ~ 2 
massive galaxies. Objects with BzK > —0.2 are classified as star-forming galaxies (sBzKs). Sources with 
BzK < —0.2 and (z — K)ab > 2.5 are defined as passive galaxies (pBzKs). 

Following the BzK technique in Daddi et al. (2004), many groups selected large samples of BzKs 
(include sBzKs and pBzKs) from different surveys (Kong et al. 2006; Lane et al. 2007; Blanc et al. 2008; 
Hartley et al. 2008; Hayashi et al. 2009; Cassata et al. 2010; McCracken et al. 2010; Onodera et al. 2010; 
Fang et al. 2012; Ryan et al. 2012). Furthermore, they also investigated the physical properties of these 
galaxies, e.g., surface density, stellar mass, SFR, near-infrared (NIR) spectroscopy, morphology, clustering, 
and size. So far the largest sample of BzKs is from McCracken et al. (2010), they found that the clustering 
of BzKs is much stronger than that of full Ff-limited (F^ab < 22) samples of field galaxies. Moreover, the 
comoving correlation length of pBzKs (rg ~ 7 h“^Mpc) is larger than that of sBzKs (rg ~ 5 h~^Mpc). 

For galaxies at z ^ 2, Hubble Space Telescope (HST)AVide Field Camera 3 (WFC3) NIR imaging can 
provide crucial clues to the rest-frame optical morphologies. At such redshift HST/WFC3 NIR bands move 
beyond the Balmer break (Arest ^ 4000 A) to the redder wavelengths and thus probe the light from the 
dominant stellar population of galaxy. This will enable us to study the rest-frame optical morphologies and 
structures of BzKs at z 2. By using HSTAVFC3 NIR images in the Hubble Ultra Deep Field (HUDF), 
Cassata et al. (2010) reported the structural features of 6 pBzKs, these galaxies in appearance have a rel¬ 
atively regular morphology and smaller size than local ellipticals of analogous stellar masses. The similar 
results were also found in Ryan et al. (2012), their pBzHs sample includes 30 passive galaxies at z ~ 2 (use 
the FF-band filter replaces the FL-band filter). Fang et al. (2012) measured nonparametric morphological 
parameters of BzKs (50 pBzKs and 173 sBzKs) at z ^ 2 in the Extended Groth Strip (EGS) field, and 
found that BzKs have both early and late types. 

This paper will utilize HST/WEC3 E160W images to investigate the structural properties of BzKs. 
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and for the first time we study the size evolution of sBzKs. We introduce the multi-band observations and 
data reduction of the COSMOS field in Section 2. Section 3 describes the selection, redshift distribution 
and number counts of BzKs. We show the SFR — M* correlation of BzKs in Section 4. We present the 
structural and morphological properties of BzKs in Section 5 and 6, and summarize our results in Section 7. 
Throughout this paper, we adopt a standard cosmology Hq = 70 km s“^ Mpc“^, = 0.7, and JIm = 0.3. 
All magnitudes use the AB system unless otherwise noted. 

2 OBSERVATIONS AND DATA 

The Cosmic Evolution Survey (COSMOS) is targeted on a special area of the sky, that has been observed 
with some of the world’s most powerful telescopes on the ground and in space, in wavelength range from X- 
rays through ultraviolet and visible light, down to infrared and radio waves (Scoville et al. 2007). In addition, 
it also includes optical/infrared spectroscopy survey using the Keck DEIMOS and LRIS, Magellan IMACS, 
and VET VIMOS spectrographs. More details of the observation and data reduction in the COSMOS field 
can be found in McCracken et al. (2012) and Muzzin et al. (2013). The multi-band photometry data we 
use in our work is from the AT-limited (K < 23.4) catalog of the COSMOS/UltraVISTA field provided by 
Muzzin et al. (2013), which is produced based on the NIR data from the Ultra VISTA DRl (McCracken 
et al. 2012). Moreover, the derived physical parameters we use in our study also come from the catalog of 
Muzzin et al. (2013), such as stellar mass (M*), SER (SFRuv.uncorr + SFRir), and photometric redshift 
(zp, if there is no spectroscopic redshift available). 

In this paper, we will utilize the latest released data of HSTAVFC3 F160W high-resolution images 
in the CANDELS*-COSMOS field to investigate the morphological features of BzKs in our sample. The 
CANDELS/wide COSMOS survey covers a total of ^210 arcmin^ at J— and 77—band. The 5 cr limiting 
magnitude is 26.9 in the E160W filter. HSTAVEC3 E160W images were drizzled to 0".06 pixeP^. Eurther 
details are in Grogin et al. (2011) for the survey and observational design, and Koekemoer et al. (2011) for 
the data products. 

3 SELECTION, REDSHIFT DISTRIBUTION AND NUMBER COUNTS OF BZKS 

To construct a sample of galaxies at z 2, we use the BzK color criteria from Daddi et al. (2004). 
The optical B- and z-band data is taken with Subaru/SuprimeCam (Bj, z’*'), while the K data is from 
the VISTA/VIRCAM. In addition, we also apply the color correction used by McCracken et al. (2010) to 
Bj — z~^ to keep unanimous with the BzK selection technique. As shown in Figure 1, blue dots represent 
14550 sBzKs with BzK > —0.2 (solid line) and K < 22.5, and red dots correspond to 1763 pBzKs with 
BzK < —0.2 and z — K > 2.5 (dot-dashed line). Objects with z — K < 0.3(77 — z) — 0.5 (dashed line), 
are classified as stars. 

Figure 2 shows the redshift distribution for BzKs in the COSMOS field. For a sample of galaxies with 
K < 22.5, in Figure 2(a), we find that the BzK color technique successfully selects more than 80% of 
galaxies at redshift 1.5 < z < 2.7 (this fraction goes to more than 90% at 1.6 < z < 2.6), indicating the 
BzK criteria is quite effective in selecting galaxies at z ^ 2. For the massive galaxies (M* > IO^^Mq) 
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Fig.l BzK two-color diagram for all objects in the COSMOS field. Galaxies with BzK > —0.2 
(solid line) and K < 22.5 are selected as sBzKs. Sources with BzK < —0.2 and z — K > 
2.5 (dot-dashed line) are defined as pBzKs. The color criterion of star and galaxy separation is 
z — K = 0.3{B — z) — 0.5 (dashed line). 


at 1.6 < z < 2.6, the percentage of objects selected as BzKs is > 90% (the dotted red lines in Figure 
2(a)). The mean redshifts of sBzKs (Figure 2(b)) and pBzKs (Figure 2(c)) are 1.75 ± 0.48 and 1.69 ± 0.33, 
respectively. In the meantime, we use the extrapolation of the red-blue separation method of Bell et al. 
(2004) to separate the red sequence and blue cloud in the rest-frame U — V vs. My diagram. About 91% 
of pBzKs can be roughly divided into red sequence, and most of sBzKs (84%) are located in the blue cloud 
region. 

We calculated differential AT-band number counts for all sBzKs and pBzKs in our sample (see Figure 
3). For comparison with the results from previous works, we also plotted their data in Figure 3 (Deep3a- 
F and Daddi-F of Kong et al. 2006; Lane et al. 2007; Blanc et al. 2008; Hartley et al. 2008; McCracken 
et al. 2010; Fang et al. 2012). The dot-dashed lines represent the counts of quiescent galaxies (QGs) and 
star-forming galaxies (SFGs) from the semi-analytic model (Kitzbichler & White 2007). In general, our 
counts agree with the results from the literature. Owing to the existence of photometric offsets and cosmic 
variance, there is a discrepancy among different works for sBzKs and pBzKs counts. Combining with the 
data from the literature, we confirm that the number counts of pBzKs has a break at K ^ 21.0, a possible 
explanation for this lies in the small redshift range of pBzKs (Kong et al. 2006). The redshift distribution of 
our pBzKs sample also supports the finding, compared to sBzKs. Compared with the observed results, the 
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Fig. 2 (a) Fraction of BzKs (include sBzKs and pBzKs) in the total galaxy sample of K < 22.5. 
The dotted red lines show the fraction of massive galaxies (M* > 1O^°M0) selected as BzKs. (b) 
Redshift distribution for sBzKs in our sample, (c) Redshift distribution for pBzKs in our sample. 
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Fig. 3 Differential number counts of sBzKs and pBzKs in the COSMOS/UltraVISTA held. The 
results from the literature and the model are also shown in this diagram. Deep3a-F and Daddi-F 
from Kong et al. (2006). 


4 THE STELLAR MASS-5Fi? CORRELATION OF BZKS 

For local star-forming galaxies, Brinchmann et al. (2004) found that there is a tight correlation between 
M* and SFR (SFRocM“), and called it main sequence (MS). At redshift 0.5 < z < 3, the MS is also 
conhrmed (Daddi et al. 2007; Elbaz et al. 2007; Rodighiero et al. 2011; Fang et al. 2012), but the slope 
(a) ranges from 0.6 to 1.0 (relying on different sample and the approaches for calculating M* and SFR). In 
Figure 4, we show the relation of M* and SFR of BzKs in the COSMOS held. For sBzKs, a best-ht slope 
a = 0.67 ± 0.06 (blue line) is found, in agreement with those provided by Daddi et al. (2007) (a ^ 0.9, 
gray line) and Rodighiero et al. (2011) {a ^ 0.79, cyan line). The discrepancy for different slopes is due to 
the different methods in deriving SFR (Rodighiero et al. 2014). On the other hand, we hnd that the SFRs 
and stellar masses of pBzKs also show a correlation, but with lower SFRs compared to sBzKs for a given 
stellar mass. Gray squares in Figure 4 represent sBzKs in Pannella et al. (2009). From Karim et al. (2011), 
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Fig. 4 Relationship of stellar mass vs. SFR for BzKs in the COSMOS held. Solid squares and 
triangles represent the star-forming galaxies from Pannella et al. (2009) and Karim et al. (2011), 
respectively. Gray and cyan lines correspond to the MS from Daddi et al. (2007) and Rodighiero 
et al. (2011), respectively. 

5 STRUCTURES OF BZKS 

In order to analyze the structural properties of BzKs in the CANDELS-COSMOS held, we employ the 
latest catalog” (version 1.0) provided by van der Wei et al. (2012). Galaxy sizes (re) are measured from 
the HSTAVFC3 F160W imaging. Within a matched radius of 0".5 , we obtain the structural parameters of 
52 pBzKs and 378 sBzKs from van der Wei et al. (2012). The M* — re relations are shown in Figure 5 of 
pBzKs and sBzKs, respectively. Shen et al. (2003) provided the relations with Itr dispersion for local late- 
and early-type galaxies (FTGs and ETGs) as plotted in this hgure. And for comparison, other data from the 
literature are also shown in Figure 5(b) (Cassata et al.2010; Gobat et al. 2012; Ryan et al. 2012; Szomoru et 
al. 2012). From Figure 5, we hnd that the sizes of pBzKs and massive sBzKs at z ^ 2 are smaller than their 
local counterparts at a hxed stellar mass. Moreover, we also see a diversity of structural properties among 
BzKs, some sources are similar to local galaxies, but there is also the existence of massive compact BzKs, 
compared to present-day counterparts. Generally, pBzKs have elliptical-like compact structures with low 
Te, while sBzKs are relatively extend and irregular, with higher rg. 

To further investigate the size evolution with redshift for our BzKs sample at z ~ 2. We show the sizes 
for pBzKs (1.85 ± 1.09 kpc) and sBzKs (2.63 ± 1.36 kpc) in Figure 6, respectively. The effective radii 
of QGs and SFGs from the literature are also plotted in this hgure (Shen et al. 2003; Cassata et al.2010; 
Gobat et al. 2012; Ryan et al. 2012; Szomoru et al. 2012; Fan et al. 2013; Patel et al. 2013; Morishita et al. 
2014). As shown in Figure 6, we hnd that the sizes of pBzKs are a factor of 2 — 3 smaller than those of 
local counterparts. For sBzKs, the mean size is one to two times smaller than those of typical local FTGs 
with comparable mass. Combined with the data points of Patel et al. (2013) and Morishita et al. (2014), the 
difference of sizes for pBzKs and sBzKs indicates that the two classes have different evolution processes 
and assembly histories, such as minor mergers with a low increase in galaxy’s mass and secular evolution 
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Fig. 5 Distribution of M* vs. for BzKs ((a) sBzKs and (b) pBzKs) in the CANDELS- 
COSMOS field. The results from the literature are also shown in Figure 5(b) (Cassata et al.2010; 
Gobat et al. 2012; Ryan et al. 2012; Szomoru et al. 2012). 



Fig. 6 Evolution of the effective radius (re) with redshift for pBzKs and sBzKs in our sample. 
The effective radii of QGs and SFGs from the literature are also plotted in this figure. The mean 
size of pBzKs is 1.85 ± 1.09 kpc (red solid circle), while sBzKs is 2.63 ± 1.36 kpc (blue solid 
circle). Left: Green and cyan lines correspond to re oc (1 + (Patel et al. 2013) and 

re oc (1 + (Morishita et al. 2014), respectively. Right: Green and cyan lines represent 

re oc (1 + 2:)“0'63 (Patel et al. 2013) and re oc (1 + (Morishita et al. 2014), respectively. 


without mergers (or monolithic collapse mode). For the size evolution of pBzKs, our observations support 
that the predictions from minor mergers. As regards the size growth of sBzKs, a possible explanation lies 
in internal evolution through some quenched mechanisms (AGN feedback or stellar winds). 

6 MORPHOLOGIES OF BZKS 


Using the WFC3 on board the HST, CANDELS provided high resolution F160W imaging data 
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Fig. 7 Distribution of BzKs in the M 20 vs. Gini coefficient plane. The red and blue circles 
represent pBzKs and sBzKs, respectively. 

of z ~ 2 BzKs. To clearly analyze their structural features, we calculated nonparametric morphological 
parameters of galaxy, such as Gini coefficient (G; the relative distribution of the galaxy pixel flux values) 
and M 20 (the second-order moment of the brightest 20% of the galaxy’s flux) (Abraham et al. 1996; Lotz 
et al. 2004). 

In Figure 7, the red and blue circles represent pBzKs and sBzKs, respectively. Meanwhile, based on the 
values of G and M 20 of 52 pBzKs and 378 sBzKs, their stamp images (3" x 3") are also plotted in this fig¬ 
ure. As shown in Figure 7, we find that pBzKs in appearance have regular and compact (like spheroid), and 
that show low M 20 and high G in rest-frame optical morphology. For sBzKs, there is a wide range of mor¬ 
phological diversities, including clumpy, irregular, extended, and early-type spiral-like morphologies, but 
most of them show diffuse structures, with high M 20 and low G. That indicates the more concentrated and 
symmetric spatial extent of stellar population distribution in pBzKs than sBzKs. Furthermore, we derived 
the mean values of G and M 20 for pBzKs and sBzKs, corresponding to (0.63, —1.70) and (0.51, —1.49), 
respectively. Our findings further imply that passive galaxies and star-forming galaxies have different evo¬ 
lution modes and mass assembly histories. 

7 SUMMARY 


Based on a SziT-selected technique, we present 14550 star-forming galaxies (sBzKs) and 1763 passive 
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field. Utilizing data from HST WFC3/F160W imaging in the CANDELS-EGS field, we investigate the 
morphological and structural diversities of these galaxies. Our findings are as follows: 

(1) We find that the Szif color technique successfully selects > 80% galaxies at redshift 1.4 < z < 2.7 
(this fraction > 90% at 1.6 < z < 2.6), indicating the BzK criteria is a quite effective galaxy selection 
method at z ~ 2. Moreover, for the massive galaxies (M* > 1O^°M0) at 1.6 < z < 2.6, the percentage of 
sources selected as BzKs is > 90%. 

(2) The differential number counts of sBzKs and pBzKs agree with the results from the literature. 
Compared with the observed results, models for galaxy formation and evolution provide too few (many) 
galaxies at high (low)-mass end. 

(3) We find that the star formation rate (SER) and stellar mass of sBzKs follow the relation of main 
sequence (SERocM° ®^^° °®). pBzks also shows correlation between SERs and stellar masses but with too 
low SERs for pBzKs. Moreover, about 91% of pBzks can be roughly divided into red sequence. 

(4) We find that the sizes of pBzKs (1.85 ± 1.09 kpc) and sBzKs (2.63 ± 1.36 kpc) at z ~ 2 are smaller 
than their local counterparts at a fixed stellar mass. Moreover, we also see a diversity of structural properties 
among BzKs, some sources are similar to local galaxies, but there is also the existence of massive compact 
BzKs, compared to present-day counterparts. The sizes of sBzKs are larger than pBzKs in general, even in 
high-mass systems, but some have very compact structures, with < Ikpc. 

(5) We find a wide range of morphological diversities for sBzKs, from extended or diffuse to early-type 
spiral-like structures, while pBzKs are relatively regular and compact (like spheroid). Moreover, we calcu¬ 
late the mean values of G and M 20 forpBzKs and sBzKs, corresponding to (0.63, —1.70) and (0.51, —1.49), 
respectively. The sBzKs show high M 20 and low G, which indicates less concentrated and symmetric spa¬ 
tial distribution of the stellar mass of sBzKs at z ~ 2, comparing to pBzKs. Our findings imply that the two 
classes have different evolution modes and mass assembly histories. 
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